Reducing number fluctuations of ultracold atomic gases via dispersive interrogation 
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We have used dispersive laser probing to follow the central density evolution of a trapped atomic cloud during 
forced evaporative cooling. This was achieved in a heterodyne detection scheme. We propose to use this as a 
nondestructive precursor measurement for predicting the atom number subsequent to evaporation and provide 
a simple experimental demonstration of the principle leading to a conditional reduction of classical number 
fluctuations. 
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Ultracold gases, notably Bose-Einstein condensates 
(BECs), have received considerable interest in recent years as 
a source for metrological applications such as atom interfer- 
ometers [1]. In this context, optical dispersive measurements 
have been put forward as a method for reducing condensate 
number fluctuations in feedback protocols [2, 3] and as a 
method for quantum state engineering of atomic ensembles 
[4]. The very act of measuring the phase shift of light 
passing through a cloud of atoms can define a collective 
pseudo-spin component beyond the standard quantum limit 
while retaining a significant atomic coherence. Indeed, 
several recent experiments have reported metrologically 
relevant measurement induced spin squeezing (see, e.g., [5] 
for a review). 

Outside the realm of measurements and stabilization at 
the quantum level, dispersive detection techniques have also 
found interesting applications. Phase contrast imaging has 
for instance been applied to film the formation of magnetic 
domains in spinor BECs [6]. Furthermore, dispersive prob- 
ing was used to follow breathing [7] and center-of-mass os- 
cillations [8] for trapped atomic clouds, real time recording 
of Rabi oscillations [9-11], and the detection and correction 
of random rotations of a collective pseudo-spin [12]. In the 
context of the present paper, it is of notable interest that dis- 
persive probing holds the potential to increase the duty cycle 
of a repeated experiment, which has direct metrological im- 
pact. For example, "nondestructive" measurements have been 
employed to reduce the Dick noise of an optical lattice clock 
by minimizing the dead time of the clock cycle [13]. 

In this Brief Report, we demonstrate the reduction of (clas- 
sical) number fluctuations of a magnetically trapped atomic 
gas after forced rf evaporative cooling to the ultracold domain. 
By using off-resonant probe laser light to measure the peak 
atomic column density near the beginning of the evaporation 
sequence we can predict with high likelihood the number of 
atoms at the end of the sequence. Conditional on the outcome 
of the dispersive precursor measurement, it becomes possible 
to interrupt the experimental cycle at an early stage in the case 
of an unfavorable starting point, thereby saving precious ex- 
perimental uptime. 

Our experimental apparatus and its operation is standard 
[14]. Initially, a 87 Rb magneto-optical trap (MOT) is loaded, 



followed by a compressed MOT and molasses stage. The 
atoms are optically pumped to the \F = 2,mp = 2) hyper- 
fine substate and transferred into a magnetic quadrupole trap 
mounted on a motorized translation stage. This stage shut- 
tles atoms to an ultra-high vacuum science cell where they 
are loaded into a "clip-style" Ioffe-Pritchard trap [15], char- 
acterized by axial and radial trapping frequencies of co fl = 
2% x 16 Hz and co r = 2% x 153 Hz, respectively. Finally, 
forced rf evaporation [16] cools our sample to a temperature 
of ~ 3 juK. 

In many applications using ultracold atomic clouds, a stable 
final atom number from run to run is desirable. Figure 1 shows 
a typical distribution in final atom numbers achieved in our 
experimental setup when continuously running at a 140 s duty 
cycle which allows for evaporation, the magnetic field coils to 
cool down, and the MOT to build up to a saturated level. The 
overall number distribution has a mean value of 1.15 x 10 6 
atoms and a standard deviation of 1.4 x 10 5 . Obviously, the 
spread in numbers could be reduced by using a protocol that 
would discriminate against initial samples likely to result in a 
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FIG. 1. (Color online) Histogram for number of atoms as observed 
in 23 consecutive experimental runs (light orange). Superimposed 
is the resulting distribution after filtering according to a dispersive 
precursor measurement (dark red). 
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FIG. 2. (Color online) Concept of measurement scheme, (a) Dispersive probing: during the evaporative cooling stage the atoms are probed 
using an off-resonant laser light frequency triplet (inset). The phase-shift and Faraday rotation of the red sideband (r) upon passage through 
the sample is represented by phasor diagrams for polarization states parallel to the carrier (C) in a frame co-rotating with C. The photodetector 
experiences a beat note at frequency £1 as a result of the rotation and shortening of r. (b) The time evolution in the dispersive signal for a low 
(left) and a high (right) number sample, (c) At the end of the evaporation atoms are probed using conventional destructive absorption imaging, 
(d) Absorption images acquired at 4.5 ms time-of-flight, corresponding to the low (left) and high (right) number (N) samples followed in (b), 
respectively. 



significant deviation from some target value. Figure 2 illus- 
trates how this can be achieved via dispersive interrogation of 
the atoms using the technique of fm spectroscopy [17, 18]. 

Our dispersive probing scheme is similar to that of 
Ref. [11]. The probe beam is generated with a diode laser and 
blue-detuned from the 87 Rb D2 line F = 2 F' = 3 transi- 
tion by 707 MHz. An electro-optic phase modulator imprints 
a red (r) and a blue (b) sideband shifted by £1 = 500 MHz from 
the carrier (C). This frequency triplet then passes through an 
acousto-optic modulator, used for pulsing the probe beam, 
which shifts the frequencies down by 80 MHz. The prob- 
ing sideband (i.e. the red sideband) is hence blue-detuned by 
A = 127 MHz from the F = 2 — >> F' = 3 transition [see inset 
panel of Fig. 2(a)]; it contains 7% of the total beam power. 
The dispersive probe beam propagates along the axial direc- 
tion (z-axis) of our atomic sample and is centered on it radi- 
ally, as illustrated in Fig. 2(a). The probe beam has a waist 
of wo = 30pm and is linearly polarized along the x-axis, per- 
pendicular to the 2.1 G bias field (pointing in the z-direction) 
at the magnetic trap center. On passage through the atomic 
sample the probing sideband incurs a dispersive phase shift 
and its polarization will rotate as a result of the Faraday ef- 
fect [19], both dependent upon the number of atoms in the 
\F = 2,m/r = 2) state inside the probe volume (the light-atom 
interaction for the carrier and upper sideband, being detuned 
by 627 MHz and 1127 MHz, respectively, from the 2 3' 
resonance, is negligible in comparison). Consequently, an op- 
tical beat signal is generated at the modulation frequency of 
£1 = 500 MHz, which we detect with a high speed biased pho- 
todetector. The resulting signal is amplified, demodulated at 



500 MHz, lowpass filtered at 100 kHz, and sampled using a 
digitizer in 200 jus time segments around the probe pulses, 
having a duration of 20 jus. 

To characterize our scheme, we probe our atomic sam- 
ple, while evaporating, at a repetition rate of 1 Hz, using 50 
light pulses. The total optical power of the probe beam is 
10 juW and the probing sideband thus contains ~ 5.5 x 10 7 
off-resonant photons per pulse. The evaporation sequence 
consists of a series of four rf sweeps (linear in time) using 
(relative) drive powers as indicated in Fig. 2(b). At the end 
of evaporation the ultracold cloud is probed using standard 
absorption imaging [20] with a resonant probe beam propa- 
gating along the j-axis as illustrated in 2(c). From this we 
determine the total number of atoms. We terminate the disper- 
sive probe pulse train after the first 50 seconds of evaporation 
as the fragility of the atomic cloud increases substantially in 
the final two evaporation ramps and destructivity becomes sig- 
nificant; when probing for the full 66 s of the evaporation se- 
quence, the sample suffers a ~ 60% loss in final atom number. 
For our laser parameters and at the center of the beam (i.e, at 
the peak intensity) the Kramers-Heisenberg formula [21] stip- 
ulates a ~ 35% loss of atoms from the magnetically trapped 
\F = 2, mp = 2) state during a probe pulse due to spontaneous, 
inelastic Raman scattering of photons from the probing side 
band as well as the carrier light [22]. The remaining, trapped 
atoms (again, taken at the peak intensity of the probe light) 
would have scattered (elastic Rayleigh process) 5 photons on 
average - each heating an atom by twice a photon recoil energy 
C^rec/kB ^180 nK) [23]. Hence a probe pulse of the present 
parameters will significantly affect an ultracold sample with a 
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radial extent r < wo. At the time of our final probe pulse the 
atomic sample has a temperature of > 30 juK and a radius of 
r> 80 jum , i.e., still reasonably robust against probing [24]. 
The associated atomic loss and heating from probing is com- 
pletely negligible at the beginning of the evaporation, where 
we are dealing with millimeter- sized samples (r ^> w Q ) con- 
taining billions of atoms at temperatures much higher than the 
recoil temperature for spontaneous emission. 

A processed data set {di} for the series of pulses i= 1...50 
is obtained by integrating over the width of each pulse and 
subtracting from this a background value. Figure 2(b) shows 
examples of the sequential development in di for two experi- 
mental runs corresponding to a high and a low number atomic 
sample. Figure 2(d) shows the corresponding absorption im- 
ages as acquired at the end of evaporation, {di} is a measure 
of the central column density of the sample along the z-axis 
as long as wo is smaller than the sample radius. The observed 
behavior in Fig. 2(b) is explained well by a simple theoretical 
model based on the spatial sample-probe overlap integral [25]. 
Initially, di is observed to increase as a result of evaporation 
and an increase in optical depth (OD), but eventually levels 
off and decreases as the radial extent of the atomic sample 
becomes comparable to and smaller than wo. 

For our characterization, we acquired data in 53 consecu- 
tive experimental runs (designated by k) where the atom num- 
bers were deliberately varied to yield a significant dynami- 
cal signal range. This was achieved by adjusting the optical 
power of the MOT laser beams, resulting in a variation in the 
magnetic trap loading. In Fig. 3 we plot the final atom number 
a(k) as inferred from absorption imaging versus the dispersive 
signal strength for the 10 th pulse d\o(k). Previous investiga- 
tions have shown a close to linear relation beween MOT levels 
and final atom number for a fixed evaporation sequence [26] . 
Similarly, we find that the atom number at the end of the evap- 
oration sequence scales linearly with the dispersive signal 10 s 
into the evaporation sequence. As a measure of the correlation 
we compute the data sample Pearson correlation coefficient 

HLMk)-{°mW-{di)} (1) 

\JlLx W) - {a^ZU Wk) - (di)} 2 ' 

Figure 4 shows the progression of r ; with dispersive probe 
pulse number i. As would be expected, the correlation coeffi- 
cient steadily increases with i as the time separation between 
di{k) and a(k) decreases. However, even for low values of i 
we note a significant correlation. For example, for the data of 
Fig. 3 (based on i= 10) we obtain no = 0.90 

The high degree of correlation immediately implies that a 
dispersive precursor measurement can be exploited to reduce 
the number fluctuations in the final sample. Whilst collecting 
samples for the number distribution in Fig. 1 we simultane- 
ously acquired the complementary dispersive data. We now 
choose the magnitude of the dispersive signal at the 10 pulse 
mark as a precursor measurement for predicting the final atom 
number and we take 1.25 x 10 6 atoms as the target value for 
our filtering process. The target number value in turn defines 
a corresponding dispersive signal target [defined by a linear fit 
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FIG. 3. (Color online) Number of atoms a(k) as measured in ab- 
sorption imaging at the end of the evaporation sequence versus the 
dispersive interrogation pulse dio(k) (•); red line is a least squares 
fit to the data, which is used to define target from a target atom num- 
ber. The dataset corresponding to Fig. 1 ( ) is decimated (■) by 
excluding points outside a 10% band around target- 



to {dio(k),a(k)} — See Fig. 3]. By excluding experimental 
realizations with a precursor value outside a 10% band around 
^target we get the filtered data shown in dark red in Fig 1 and 
the inset of Fig. 3, superimposed on the raw data. The filtered 
number distribution has a mean of 1.27 x 10 6 and has a width 
characterized by a standard deviation of 6.1 x 10 4 which is a 
reduction of factor of > 2 in comparison to the original unfil- 
tered distribution. 

In conclusion, we have demonstrated that the information 
gained by probing a trapped atomic sample at the beginning 
of the evaporation sequence using a reasonably simple non- 
destructive, dispersive heterodyne scheme displays a signif- 
icant correlation with the final resulting atom number. This 
may be used in a preselecting protocol to achieve consistent, 
high number samples. While the precursor interrogation pulse 
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FIG. 4. Correlation coefficient r ; versus pulse number i. 
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only samples a fraction of the cloud via its column density at 
a particular point it provides a quite robust quality measure. 
Such complimentary diagnostics tools may prove valuable as 
cold atoms experiments move towards applications outside the 
laboratory [27]. We also note that nondestructive probing of 
an atomic sample during evaporative cooling may open up 
paths for more efficient adaptive optimization procedures than 
conventional and tedious algorithms using absorption imaging 
[14]. In particular, stochastic optimization of BEC production 
using genetic algorithms [28], presently based on absorption 
images as the sole input, may benefit from the additional in- 
formation gained. Along these lines, we point out that field 
programmable gate arrays would provide an ideal solution for 
acquiring and processing the dispersive signal and actuate on 



the sample in a feedback loop. 

Minimally destructive detection techniques were recently 
discussed and compared in detail by Ramanathan et al. [29] . 
One of the benefits pointed out was the possibility of normal- 
izing a resulting experimental measurement to a nondestruc- 
tive measurement on the initial sample to overcome the effects 
of number fluctuations [30]. Based on the result presented in 
this Brief Report, we argue that it should be possible to em- 
ploy our scheme to this effect by normalizing a conventional 
destructive absorption image to that of simple prior dispersive 
measurement of the peak column density of the sample. 
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